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ABSTRACT 



Method and apparatus for non-invasive testing of digital 
communications systems. Amplitude measurements are 
made for multiple frequencies of a multi-frequency com- 
munication system, converted to the time domain. An adap- 
tive filter output is matched to the time domain representa- 
tion to characterize the channel. Impedance mismatches may 
be precisely located using this technique. An error signal 
representing a difference between a signal transmitted 
through the channel and a received signal is estimated and 
analyzed. The error signal is separated into components 
corresponding to contributions by wide band noise, residual 
phase modulation, and residual amplitude modulation. Iden- 
tification and removal of narrow-band interferes may occur 
prior to this separation. Bit error rate and system margin 
computations employ a Monte Carlo simulation of the 
various error sources. This provides a well refined estimate 
of bit error rate and system margin. 

26 Claims, 16 Drawing Sheets 
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NON-INVASIVE DIGITAL CABLE TEST estimated and analyzed. The error signal is separated into 

SYSTEM components corresponding to contributions by wide band 

noise, residual phase modulation, and residual amplitude 

STATEMENT OF RELATED APPLICATIONS modulation. Identification and removal of narrow-band 

_ . .. . irn^ xt 5 interferers may occur prior to this separation. Bit error rate 

This application is a division of U.S. app. Ser. No. , . r . r , w*^, 

„ n -J/-. . f ~ n . orv7 tt c n * m in* an d system margin computations employ a Monte Carlo 

08/791386, filed Jan. 30, 1997, now U.S. Pat. No. 6,061,393 • , %• *u * • i * *u m.- 

, . , 7 . . . . r TT ^ ' ' simulation that simulates the vanous error sources. This 

which is in turn a continuation-in-part of U.S. patent app. „ c , t cu .. . , 

" KT £1 . 0 ^ ¥1 v „ . /J provides a well refined estimate of bit error rate and system 

Ser. No. 08/532,017, filed Sep. 21, 1995, now U.S. Pat. No. margm 

5,751,766 which is in turn a continuation-in part of U.S. . " , ., . . . *■ i ■ 

patent application Ser. No. 08/430,909, filed Apr. 27, 1995, 10 , In accordance with a third aspect of the invention^ a 

now abandoned. The contents of both parent applications are cbi f { variance measurement based on an error sigpal is 



herein incorporated by reference for all purposes. 



used to estimate bit error rate. In accordance with a fourth 
aspect of the invention, an error signal is displayed in the 

BACKGROUND OF THE INVENTION frequency domain relative to a measured power of a 

15 received signal. 

The present invention relates to measuring signal degra- ^ invention ^ ^ better by refererice t0 

datum in commumcation systems. In particular, the inven- ^ following detailed description in connection with the 

tion relates to methods and apparatus for non-invasively accompaaying favnnss. 
measuring performance of a digital broadcast system. 

Digital communication technologies offer numerous 20 BRIEF DESCRIPTION OF THE DRAWINGS 

advantages over their analog predecessors. Nearly error-free FIG. 1 depicts a simplified representation of a represen- 

transmission is guaranteed if a threshold signal-to-noise is tative digital communication system suitable for testing in 

available. Channel distortions may be corrected using adap- accordance with the present invention, 

live equalizers. Coding techniques may be advantageously FIG. 2Ais a block diagram depicting one implementation 

employed both to overcome channel-related signal impair- 25 of a portion 0 f a digital receiver. 

ments and to minimize usage of bandwidth. nG 2fi whkh depicts a sigQal constellation for 

Nonetheless, digital commumcation system operation is a representative digital modulation scheme, 64 point QAM. 

not fault-free. The transmission of information over a physi- FIG 2C depicts an alternative implementation of a portion 

cal medium ultimately requires along components such as Q ^ a digital receiver 

RF mixers, amplifiers, oscillators, etc. that are prone to 30 „ T ^, - . . " , . f 

. v * + j 1 *r>. j - „+ u ~~ FIG. 3 is a simplified representation of a diagnostic 

misalignment, temperature caused drift, and vanous other . . . K , y 

j rc -i rJ[ , . t ir 1 - processor/momtor in accordance with one embodiment of 

modes ot railure. The physical medium itself also introduces , 

4 . iU / £ .j j • j a the present mvention. 

impairments in the iorm 01 added noise and reflections. r 

Digital receivers can inherently correct for errors up to a FIG 4A depicts hardware for implementing a communi- 

certain threshold, but if that threshold is exceeded, commu- * catlon diagnostic system in accordance with a preferred 

nication is severely impaired. embodiment of the present mvention. 

Digital communication techniques found their earliest JJ^^**** * SyS,6m bl °° k ° f 3 

application in the context of point- to -point links as would be * . 

used by the military or commercial common carriers. In FIG. 5 , A 15 a flowchart describing the steps of analyzing 

these applications, the link may be tested by interrupting 40 a channel response in accordance with one embodiment of 

normal service and transmitting test signals. For example, if ™ P resent m vcntl011 - 

it is desired to test the link margin of a link, a test signal FIG 5B shows a svstem for implicitly calculating the 

could be transmitted at progressively reduced power levels pseudo-inverse of an equalizer response in accordance with 

until the signal can no longer be accurately received. one embodiment of the present invention. 

Increasingly, digital techniques are being applied to 45 5C is a flowchart describing steps of analyzing a 

broadcasting. One important example is broadcasting of channel response by making measurements at multiple fre- 

digital video signals to numerous subscribers over coaxial quenacs in accordance with one embodiment of the present 

cable, or over the air. Commercial digital broadcast systems invention. 

require high reliability yet these digital broadcast systems 5Q FIG - 5D depicts a typical structure of a channel in the 

cannot be tested as easily as the point-to-point links. Inter- context of an HFC (hybrid fiber cable) system, 

rupting normal service for transmission and measurement of FIG. 6A is a block diagram describing various modes of 

test signals is not a commercially viable option. taking advantage of error signal 217 which represents 

ingress interference in communications channel 104. 

SUMMARY OF THE INVENTION ^ FIG. 6B depicts a communications channel with an addi- 

The present invention provides method and apparatus for ti °nal ingress signal being injected by a second digital 

non-invasive testing of digital commumcation test systems. transmitter. 

In accordance with a first aspect of the invention, amplitude FIGS. 7A-7B are diagrams illustrating how system mar- 
measurements are made for multiple frequencies of a multi- gin may be tested non-invasively in accordance with one 
frequency communication system, and converted to the time so embodiment of the present invention, 
domain. An adaptive filter then seeks to match its output to FIG. 7C is a diagram illustrating the effect of broadband 
the time domain representation, thereby characterizing the noise on a particular constellation point as received, 
channel. Impedance mismatches may be precisely located FIG. 7D is a diagram illustrating the effect of a narrow- 
using this technique. band interferer on a particular constellation point as 

In accordance with a second aspect of the invention, an 65 received, 

error signal representing a difference between a signal FIG. 7E is a diagram illustrating the effect of residual 

transmitted through the channel and a received signal is phase noise on a particular constellation point as received. 
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FIG. 7F is a diagram illustrating the effect of residual converted signal to a discrete-time digital representation, 

amplitude modulation on a particular constellation point as The sampling rate is . synchronized to symbol boundaries 

received. within the received data stream by a baud time recovery loop 

FIG. 7G is a flowchart describing steps of analyzing an 212 that generates a sampling rate control signal 213. 

error signal to estimate system margin and obtain other 5 An adaptive equalizer 214 automatically corrects for 

system information in accordance with one embodiment of distortions in channel 104 and typically includes a digital 

the present invention. FIR fil^ 1 and/or HR filter with variable tap weights. 

FIG. 7H is a diagram depicting narrowband interference (Alternative receiver architectures employ equalizers with 

superimposed on a digital signal. fi xed . ta P weights or omit equalization entirety.) A symbol 

V™ 1, . * - • t * .10 decision stage 216 exammes the data output from adaptive 

FIG 71 is a diagram depicting the response of a notch ^ 2U ^ ^ taasa ^ vM data based orj a 

filter for removing narrowband interference in accordance maximum likelihood decision rule corresponding to the 

with one embodiment of the present invention. modulation scheme employed. A channel decoding stage 

FIG. 7J is a diagram depicting the geometry of a particular 218 removes effects of forward error correction (FEQ or 

constellation point. 15 other channel coding schemes employed. Further stages (not 

FIG. 7K is a diagram depicting a constellation point as shown) may remove source coding, differential coding, 

affected by more than one type of impairment. scrambling, interleaving, or compression applied to the 

FIGS. 8A-8C are diagrams illustrating how parameters transmitted data, 

relating to residual modulation and symbol timing variation Symbol decision stage 216 also generates an error signal 

may be measured in accordance with the invention. 20 217 corresponding to the difference between an ideal modu- 
lation signal representing the output estimates and input of 

DESCRIPTION OF SPECIFIC EMBODIMENTS symbol decision stage 216. This error signal is an input to a 

Digital Communications Background weight updating stage 220 that adjusts weights 221 of 

FIG. 1 depicts a simplified representation of a represen- adaptive equalizer 214 to minimize certain components of 

tative digital communication system 100 suitable for testing 25 the error. A carrier recovery state 222 also makes use of this 

in accordance with the present invention. Digital commu- error signal to remove residual frequency offset and fre- 

nication system 100 includes a transmitter 102, a commu- quency modulation effects induced in part by tuning/ 

nications channel 104, and a digital receiver 106. Transmit- downconversion stage 208 by controlling the internal VCO 

ter 102 receives digital data to be transmitted to receiver 106 with a desp inning control signal 223. 

and incorporates channel coding, modulation, 30 The operation of typical digital communication system 

up-conversion, amplification, and output filtering functions. 100 will be made clearer by reference to FIG. 2B which 

The present invention is not restricted in application to any depicts a so-called signal constellation for a representative 

particular channel code, digital modulation scheme, or fre- digital modulation scheme 64 point QAM. A signal constel- 

quency band. Channel 104 represents any physical medium lation is a graphical representation of the possible symbols 

for transferring a signal for transmitter 102 to receiver 106. 35 for a given modulation scheme. The horizontal and vertical 

For example, transmitter channel 104 may represent a axes correspond to the orthogonal components of the modu- 

coaxial cable delivery system and transmitter 102 may be a lation signal. Each possible symbol is represented by a point 

cable head-end. Alternatively, channel 104 may represent at the position of its associated (1, Q) coordinates. As shown 

transmission between a satellite and a ground-based in FIG. 2B, 64 point QAM is represented as an array of 64 

receiver. Furthermore, receiver 106 need not be uniquely 40 points. Since log 2 (64)-6, the choice of one particular sym- 

coupled to channel 104. Channel 104 may represent a bol for transmission during a given symbol period commu- 

broadcast medium accessible to numerous receivers. nicates 6 bits of information. 

The present invention takes advantage of various signals Channel 104 introduces noise causing the constellation 

and data generated internally to digital receive 106. FIG. 2A points to disperse from their ideal locations marked in FIG. 

is a block diagram depicting one implementation of a 45 2B. Transmitter filtering and echoing effects within channel 

portion of digital receiver 106. After filtering and 104 introduce inter-symbol interference (ISI) which must be 

downconversion, a signal received from channel 104 is removed by equalization. Also, tuner/downconversion stage 

passed through a variable gain stage 202. The amplified 208 introduces residual frequency offset which causes the 

output of gain stage 202 passes through an AGC control received signals to "spin" about the origin, 

circuit 204 which controls the gain of gain stage 202 using 50 Symbol decision stage 216 essentially receives from 

a gain control signal 205. A bandpass filter 206 selects a adaptive equalizer 214 a series of points in the coordinate 

particular range of frequencies so as to guard against later space represented by the graph of FIG. 2B and assigns to 

development of unwanted image signals by a tuning/ each point the nearest constellation point as the symbol 

downconversion stage 208. Tuning/downconversion stage received. Weight updating stage 220 adjusts equalizer 214 to 

208 selects a particular signal of interest and converts this 55 minimize the distance between the received signal and the 

signal to baseband. Typically, tuning/downconversion stage constellation point transmitted. Carrier recovery stage 222 

208 incorporates at least one voltage controlled oscillator operates to minimize the spinning of received points around 

(VCO) which is at least partially controlled by a special the origin. 

control signal developed to minimize residual frequency FIG. 2C depicts an alternative implementation of a por- 

modulation in the downconverted signal. For modulation 60 tion of digital receiver 106 as described in U.S. Pat. No. 

schemes that involve two orthogonal components such 5,495,203 issued at Harp et al., EFFICIENT QAM 

QPSK, QAM etc., tuning/downconversion stage 208 devel- EQUALIZER/DEMODULATOR WITH NON-INTEGER 

ops two outputs corresponding to each orthogonal compo- SAMPLING, assigned to the assignee of this application, 

nent. The present invention is not limited to any particular and the contents of which are herein incorporated by refer- 

modulation scheme. 65 ence for all purposes. This receiver architecture will now be 

Up until this stage, the signals are presumed to be analog described in brief. An A-D converter 224 corresponds to 

and continuous. A sampling stage 210 converts the down- sampling stage 210 of FIG. 2 A but does not operate at 
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baseband but rater accepts an input centered at an IF and From despinning control signal 223, equalizer weights 221, 
samples at a fixed rate determined by a reference oscillator and error signal 217, diagnostic processor/monitor 300 
226 in conjunction with a divide by 4 circuit 228. The fixed generates information about phase noise, residual frequency 
sampling rate and IF are chosen in relation to the spectrum modulation, frequency error, and frequency stability. From 
of the modulation signal so as to enable digital quadrature 5 sampling rate control signal 213 and equalizer weights 221, 
direct conversion to baseband by a digital QDC state 230. A diagnostic processor/monitor 300 generates information 
digital rcsampler 232 synchronizes to the received data in about transmitting baud timing. Equalizer weights 221 may 
conjunction with a symbol timing phases lock loop 234 be t0 f ne " te ^formation about the channel response 
analogous to baud timing recovery loop 212. An adaptive including the time domain locations of impedance mis- 
equalizer 236 is specially adapted to operate at a sampling 10 patches. Error signal 217 generated by symbol decision 
4 t . . . . > * *• r*L. u i * a . stage 216 can be used to recover ingress signals and obtain 
rate that is non-integer fraction of the symbol rate. A earner mformation such M dus £ r vari ^ M defined 

recovery stage 238 performs the despinning function w ^ ^ ^ ^ 

directly on the digital data output from equalizer 236 with- FIG. 4A depicts hardware for implementing a communi- 
out adjustment of oscillators as in the receiver of FIG 2A. A cations link diagnostic system 400 in accordance with a 
symbol decision stage 240 estimates the received data based 15 preferred embodiment of the present invention. Communi- 
on the output of carrier recovery stage 238 and generates cations link diagnostic system 400 includes an RF turner/ 
error signals to drive carrier recovery stage 238 and an downconverter 402, a Model 194 Wideband Signal Acqui- 
equalizer weight updating stage 242. sition and Playback Unit 404 available from the assignee of 
It can be seen from FIG. 2A and 2C, that digital receivers this application, and a workstation 406. Turner/ 
generate various signals and data that correct for impair- 20 downconverter 402 receives a signal from channel 104 and 
ments in overall communication system operation. For is particularized for the band and transmission medium of 
example, equalizers 214 and 236 correct for echoing and interest. Turner/downconverter 402 selects a signal of inter- 
frequency response distortions in channel 104. AGC loop est and downconverts it to a predetermined IF. In the 
204 inherently corrects for hum and residual amplitude preferred embodiment, the IF is 43.75 MHz. 
modulation. The operation of baud timing recovery loop 25 Wideband Signal Acquisition and Playback Unit 404 
212, symbol timing PLL 234, and carrier recovery stages emulates the operation of the receiver components depicted 
222 and 238, inherently correct for residual frequency in FIG. 2A. The IF signal is digitized at up to 250 MHz and 
modulation, phase noise, drift, and instability, whether these captured in an internal snapshot memory. The remaining 
artifacts are introduced at the transmitter or receiver end. Of functions of the receiver of FIG. 2A are performed by 
course, if these various impairments exceed a given level, 30 workstation 406. 

the receiver will no longer be able to recover the transmitted Workstation 406 includes a processing unit 408 contain- 

data at an acceptable bit error rate (BER). ing control elements. A keyboard 410 and a mouse 412 are 

In general FIGS. 2A and 2C are merely representative of coupled to processing unit 408 and enable the system user to 

possible digital receiver architectures to which the present input commands and data to workstation 406. A display 

invention is applicable. Note that some components of 35 device 413 is provided for displaying test results. Worksta- 

digital receiver 106 are not present in every implementation. tion 406 further includes a floppy disk drive 414 adapted to 

For example, some implementations use equalizers that do receive a floppy disk 416, and an internal hard drive 418. 

not require weight updating. Other implementations do not FIG. 4B depicts a system architecture block diagram of 

employ any kind of equalizer. Many links do not employ workstation 406 of FIG. 4A. The system elements are 

channel coding to improve performance and thus do not 40 interconnected by a system bus 420 and include a display 

require a channel decoding stage. adapter 422 coupled to display device 413, an I/O controller 

Overview of Diagnostic System 424, system memory 425, a central processor (CPU) 426, a 

In accordance with one embodiment of the present speaker 428, a serial port 429, keyboard 410, a hard disk 

invention, impairments are non-invasively measured and controller 432 coupled to hard disk drive 418, and a floppy 

otherwise characterized by monitoring and analyzing inter- 45 disk controller 434 coupled to floppy disk drive 434. CPU 

mediate data and signals generated internal to a digital 426 controls operation of workstation 406. Program and data 

receiver. The interpretation and availability of these signals information are stored in system memory 424, on floppy 

will of course vary depending on the internal architecture of disk 416, and on hard disk drive 418. 

the particular digital receiver employed. Analysis software stored in system memory 424, on 

FIG. 3 is a simplified representation of a diagnostic 50 floppy disk 416, and/or on hard disk drive 418 implements 

processor/monitor 300 in accordance with one embodiment the functionality of diagnostic processor/monitor 300 as well 

of the present invention. Possible inputs to diagnostic as the functionality of digital receiver 106. Of course any 

processor/monitor 300 include gain control signal 205 from tangible medium including CD-ROM, magnetic tape, flash 

AGC control circuit 204, a despinning control signal 223 memory etc. could store this software, 

from carrier recovery stage 222, a sampling rate control 55 The workstation of FIGS. 4A-4B is merely representative 

signal 213 from baud rate recovery loop 212, the current of hardware usable to implement the present invention. In an 

equalizer weights from adaptive equalizer 214, and the error alternative embodiment, the present invention is imple- 

signals and received symbols generated by symbol decision mented as a battery-operated portable system suitable for 

stage 216. The monitor component of diagnostic processor/ field work. 

monitor 300 receives all these signals and converts them to 60 In the embodiment shown in FIGS. 4A-4B, Wideband 

a form suitable for analysis. The inputs of FIG. 3 are Signal Acquisition and Playback Unit 404 captures sufficient 

described in reference to the receiver architecture of FIG. 2A signal information for a full analysis in less than one second, 

but similar inputs could be obtained from the receiver The analysis itself completes in 1-2 minutes. An optional 

architecture of FIG. 2C coprocessor card installable within workstation 406 can 

From gain control signal 205, equalizer weights 221, and 65 reduce analysis time to approximately 10 seconds, 

error signal 217, diagnostic processor/monitor 300 generates In a preferred embodiment, workstation 406 derives and 

information about residual amplitude modulation and hum. displays on display device 408 combinations of the follow- 
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ing: (1) transmission (or propagation) channel response; (2) 
carrier-to-noise ratio (CNR); (3) spectrum of the underlying 
noise and any interference signals; (4) phase noise estimates; 
(5) signal spectrum, signal polar plot, signal constellation 
and eye diagram (for analyzing intersymbol interference) 5 
and (6) other diagnostic displays. 
Analysis of Single Channel Response 

In accordance with one embodiment of the present 
invention, the response of channel 104 or a portion of 
channel 104 may be determined. From the channel response, 10 
the operator can determine undesirable multipath, cable, or 
waveguide mismatching. When channel 104 is the ambient 
atmosphere, structures such as buildings, vehicles, aircraft, 
hills, etc., can reflect the transmitted signal such that the 
transmitted signal is subject to multipath. If channel 104 is 15 
a coaxial distribution system, reflections may be caused by 
impedance mismatches resulting from loose connectors, etc. 
Receiver 106 receives multiple versions of the transmitted 
signal at slightly different times due to the differences 
between the path lengths of the reflected signals and the path 20 
length of the transmitted signal. Discovering such reflections 
is an important aspect of diagnosing a link. 

The channel response analysis technique of the present 
invention will be made clear by a simplified mathematical 
characterization of digital communications system 100. The 25 
notation used herein follows a standard linear time-invariant 
model. The conversion of the results to slowly time-varying 
or somewhat non-linear systems is well-known in the art. 

Referring again to FIG. 1, transmitter 102 and associated 
electronics shape and transmits an input symbol stream S n . 30 
Channel 104 carries the output pulses corresponding to the 
symbol stream, modifying the pulses according to the 
response of the channel. Referring to FIG. 2A, the compo- 
nents of receiver 106 to the left of adaptive equalizer 214 
introduce a further response component. Hie combined 35 
response of transmitter 102, channel 104, and components of 
receiver 106 to the left of adaptive equalizer 214 can be 
denoted as H, while the response of adaptive equalizer 214 
can be denoted as W. The overall response of the system up 
to the output of adaptive equalizer 214 can then be written 40 
as 

WHS^w (1) 

where §„+e is the output of the equalizer structure and is a 
noisy estimate of the original symbol stream § M . 45 

When communications system 100 is working effectively, 
S n is a very good estimate of S„, i.e., the received symbols 
are very close to the same symbols transmitted. In this case, 
W operating on H as a matrix operator is very close to the 
identity, I, a matrix operator being defined as an operator 50 
which takes one or more vectors, matrices, or tensors as 
arguments and returns a vector, matrix, or tensor. Thus, the 
following relationships hold: 



s„~$ w (2) ss 

WHS„~S„~WH~I (3) 

It then follows that the pseudoinverse of W, W#, is very 
close to H, the overall system response. One can estimate H 
as: 60 

H-w# (4) 



FIG. 5 A is a flowchart describing the steps of analyzing 
a channel response in accordance with one embodiment of 
the present invention. At step 502, the pseudoinverse of the 65 
equalizer response is determined using the weights gener- 
ated by weight updating stage 220. 
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In one embodiment, the pseudoinverse is computed at 
step 502 explicitly, using Singular Value Decomposition as 
is well known to those of skill in the art. Any of a variety of 
well-known matrix algebraic methodologies could be used. 

FIG. SB shows a system 550 for implicitly calculating the 
pseudoinverse at step 502. An adaptive filter 552 (distinct 
from adaptive equalizer 214) receives an ideal modulation 
signal representing the symbol decisions generated by sym- 
bol decision stage 216. The output of adaptive filter 552 is 
subtracted from the input data to adaptive equalizer 214 
within a summer 556. A weight updating stage 558 then 
operates to modify the tap weights of adaptive filter 552 to 
minimize the output of summer 556, preferably using the 
well-known least mean square (LMS) procedure. Since the 
filter taps have converged, i.e., adapted to a solution such 
that the error cannot be further minimized and the error is 
acceptably small enough, the response of adaptive filter 552 
matches the overall system response. It has been found that 
the channel can be modeled with residual errors 35-40 dB 
below the level of the channel response. 

Once a pseudoinverse is computed as an estimate of the 
overall system response H, particular components of the 
system response may be isolated, at step 504. For many 
communications systems, the response of the transmitter 
102, the pre-equalizer components of receiver 106, or both 
are well known. In these cases, the known responses can be 
deconvolved from ft, leaving an estimate of the channel 
response alone. 

Furthermore, there may also be known portions of the 
response of channel 104 that may be deconvolved to focus 
on a particular portion of interest. For example, it may be 
desirable to analyze only the portion of a coaxial cable 
system lying within a subscriber's home. Once a desired 
response is obtained, it may be shown on display device 508. 

In accordance with the invention, an all-pass portion of 
the response of a channel component of interest may be 
isolated at step 506 to identify the locations of reflections. 
The time domain representation of the all-pass component 
will show the locations of individual reflections which could 
be caused by e.g., loose connectors. In accordance with the 
invention, the allpass portion of the response may be 
obtained without exact knowledge of the frequency- 
dependent amplitude response by using homomorphic signal 
processing or parametric signal processing. 

Homomorphic signal processing is a method by which a 
signal is decomposed into an all-pass component and a 
minimum-phase component. This decomposition can be 
performed using techniques well known to those skilled in 
the art, including methods described in [Oppenheim89] and 
[IEEE79], the contents of which are herein expressly incor- 
porated by references for all purposes. In the case where the 
input signal consists of a signal corrupted by micro- 
reflections, the minimum-phase portion of the signal will 
contain a representation of the spectrum of the direct arrival 
of the desired signal, while the all-pass portion will contain 
a representation of the impulse response of the channel's 
micro-reflections. 

In cases where the ideal (uncorrupted) channel response is 
known, or can be reasonably estimated, parametric signal 
processing techniques may be employed to better estimate 
the corrupted channel response. In using parametric 
approaches, a model for the corrupted channel response is 
hypothesized. This model includes several unknown 
parameters, such as time-of-arrival and amplitude of micro- 
reflections. The values for each of these parameters may 
then be determined in any number of methods well known 
to those skilled in the art, including Maximum Likelihood, 
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Maximum Entropy, Minimum Variance, etc. Details of these 
methods, as known to those of skill in the art, may be found 
throughout the literature, including [VanTrees68], 
[Haykin83], [Haykin85], [Haykin91], and [Marple87], all 
the contents of these references being herein expressly 
incorporated by reference for all purposes. 

Thus, the locations of individual reflections may be 
obtained without resort to invasive probing of the system 
while it is operating. This is an important advantage over 
prior art time domain reflectometry techniques which 
require that normal operation be interrupted. 

The resolution of the above-described system in locating 
reflection sources is limited by the effective temporal reso- 
lution of the equalizer. For a simple finite impulse response 
(FIR) structure, the effective temporal resolution is directly 
related to the tap spacing. If standard processing approaches 
are utilized, a typical T- or T/2-spaced equalizer is limited to 
a resolution of Tor T/2, respectively, where T is equal to the 
symbol period. For typical digital cable television (CATV) 
applications, a 5 MBaud, T/2-spaced equalizer would have 
a temporal resolution of 100 ns. This translates to a resolu- 
tion of about 24 meters (m) in typical coaxial cable. 

In accordance with one embodiment of the invention, by 
employing a non-integer, fractionally-spaced equalizer, the 
temporal resolution can be increased by interpreting the 
equalizer as a polyphase structure. In such systems, the 
sample rate is approximately equal to (2*N*K/L) times the 
symbol rate where 2*L>K>L, and K, L, and N are integers. 
In the case of a T/1.2-spaced equalizer, as is described in the 
Harp et al. patent, the resulting temporal resolution is T/6. 
For a 5 MBaud system, this would provide a 33 ns 
resolution, better than 8 m in the same coaxial cable where 
only 24 m was the resolution available using a T/2-spaced 
equalizer. Other spacings can provide even better resolution 
than that of a T/l .2 equalizer, though usually at the expense 
of greater computational complexity. 

With a T/1.2-spaced equalizer, one can visualize the 
structure as a linear, time-varying system comprised of five 
distinct equalizer phases. Each of these phases is used in turn 
as each new sample is received. Without loss of generality, 
the first of these phases (Phase 0) is assigned a relative delay 
of 0 and is responsible for providing an output correspond- 
ing to some symbol sent at t=0. The next phase (Phase 1) of 
the equalizer is then responsible for predicting the symbol 
sent at t=T, but using samples collected with a relative delay 
of 5T/6. Thus the taps of Phase 1 are effectively T/6 (T-5T/6) 
advanced relative to those in Phase 0. Similarly, the remain- 
ing phases are advanced by T/6 from the previous phase. 
This relationship may be summarized as in Table 1: 

TABLE 1 



Relative Delay/T 



Tbp 0 


1 


2 


3 


4 


5 


6 


7 . . . 


Phase 
















0 0/6 


5/6 


10/6 


15/6 


20/6 


25/6 


30/6 


35/6 . . . 


1 -1/6 


4/6 


9/6 


14/6 


19/6 


24/6 


29/6 


34/6 . . . 


2 -2/6 


3/6 


8/6 


13/6 


18/6 


23/6 


28/6 


33/6 . . . 


3 -3/6 


2/6 


7/6 


12/6 


17/6 


22/6 


27/6 


32/6 .. . 


4 -4/6 


1/6 


6/6 


11/6 


16/6 


21/6 


26/6 


31/6 . . . 



As shown in this table, the relative delay decreases by T/6 
from phase to phase, and increases by 5T/6 from tap to tap. 
As can also be seen from this table, by properly ordering the 
N taps of the entire equalizer structure, all relative delays 
from -4T/6 up through (N-5) T/6 in increments of T/6 are 



78,730 Bl 

10 

available. This defines the response of the equalizer with a 
resolution of T/6 (<0.17 T), three times better than that 
obtainable with T/2 equalizers and their inherent 0.5 T 
resolution. 

5 Analysis of Multi-Frequency Multi-Channel Response 

Even finer resolution in identifying reflection sources is 
obtainable by making measurements at multiple frequencies 
across the bandwidth of the entire cable system as opposed 
to monitoring a single signal. The cable system bandwidth is 
typically 200 MHz to 1 GHz. For a 500 MHz system, 
reflection sources may be localized to within 0.7 m. 

FIG. 5C is a flowchart describing steps of analyzing a 
channel response by making measurements at multiple fre- 
quencies in accordance with one embodiment of the present 
invention. At step 560, the preferred embodiment measures 

15 envelope amplitude for each signal carried by the cable 
system. Typically, a measurement is made every 6 MHz by 
tuning receiver 106 to each channel. The amplitude mea- 
surement must be made at a stage of receiver 106 prior to the 
circuitry shown in FIGS. 2A and 2C, e.g., at the input to gain 

20 stage 202. The result of step 560 is a list of frequencies and 
measured amplitudes. 

At step 562, the list of frequencies and measured ampli- 
tudes is converted into the time domain by applying an FFT 
(Fast Fourier Transform). This results in a time domain 

25 representation of the signal measured over an aggregation of 
the frequencies. At step 564, an impulse is repeatedly fed 
into an adaptive digital filter. The adaptive filter operates to 
modify its transfer function to force its output to equal the 
time domain representation developed in step 562. The 

30 adaptation method is preferably the well-known gradient 
descent algorithm. In one embodiment, the adaptive digital 
filter is an all-pole filter, and adaptation includes positioning 
only the poles of the transfer function. In another 
embodiment, the adaptive digital filter includes zeroes as 

35 well which are positioned during adaptation. 

The adaptation process forces the response of the adaptive 
digital filter to approximate the response of channel 104 as 
measured across the entire system bandwidth. This is an 
example of parametric signal processing. At step 566, the 

40 adaptive digital filter response is converted to the time 
domain. If the filter response includes both poles and zeroes, 
only the "poles" component of the response need be con- 
verted to the time domain for this purpose. The resulting 
time domain representation of the poles or "resonances 11 

45 will show pulses corresponding to the locations of imped- 
ance mismatches, as caused by, e.g., loose or corroded 
connectors. 

At step 568, if the filter response includes zeroes as well, 
the "zeroes" component of the response is isolated and 

50 displayed in the frequency domain. This shows the effects of 
shaping filters in transmitter 102, filters used at intermediate 
points in channel 104 and any other frequency response 
distortions along the way. 
Localization of Channel Impairments 

55 Generally, accurate identification of the location of an 
impedance mismatch requires isolation of a portion of the 
response of channel 104. FIG. 5D depicts a typical structure 
of channel 104 in the context of an HFC (hybrid fiber cable) 
system. A fiber-optic transmitter 570 transmits optical 

60 energy to an electro-optical node 572 Electro-optical node 
572 transmits signals through a distribution system 576 to a 
tap box 578. At tap box 578, an electrical signal is tapped off 
for transmission through drop wiring 579 to premises 580 of 
a particular subscriber. Premises 580 include televisions 582 

65 and a set-top unit 584 implementing receiver 106. 

In FIG. 5D, nodes A-E represent points where the signal 
may be tapped off and monitored to identify the channel 
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response up to that point, by any of the procedures discussed or equalizer 236) and is generally considered to be the power 

in reference to FIGS. 5A-5C. lb isolate a particular segment of the signal divided by the sum of the power of the in-phase 

of channel 104, responses are obtained at a node following error plus the power of the quadrature error. Alternatively, 

the segment and anode preceding the segment. To obtain the the cluster variance may be measured based on a modulated 

response of the optical system represented by fiber-optic 5 si g nal carrying a reencoded version of the data output from 

transmitter 570 and electro-optical node 572, the response channel decoding stage 218. The cluster variance represents 

obtained at node A is deconvolved from the response an estimate of the signal-to-distortion power rauo. Depend- 

obtained at node B. To obtain the response of distribution J5* 5? the <*f 0Ti ™ tv P e ^ verit ^ th f **l error rate 

system 576, the response obtained at node B is deconvolved £™) can be calculated as a function of the signal-to- 

/ it _ * . . 4 , „ rj, . 4U , n distortion power ratio or cluster variance. 

from the response obtained at node C. To obtain the response 10 n . r . . . . , . . . 

r , .\ „ rt . , . « j j-y . Precise calculation of the signal margin (the difference 

of drop wiring 579 the response obtained at node : C is between the signal-to-noise (SNR) ratio at the required BER 

deconvolved from the response obtained at node D. To vg me actual SNR for a ^ ^ painnent) ^ aocordance 

obtam the response of the wiring inside subscriber premises, ^ the present inveQtion is disaisse d below. However, the 

the response obtained at node D is deconvolved from the prese rjt invention also provides for a simplified calculation 

response obtained at node E. 15 0 f g£R and signal margin directly from the cluster variance 

Analysis of Error Signal based on the recognition that in most cases distortion can be 

In accordance with the invention, diagnostic processor/ considered to have an approximately Gaussian distribution 

monitor 300 receives error signal 217 corresponding to the and be independent. The Gaussian assumption holds gener- 

difference between an ideal modulation signal representing ally in digital cable systems since the equalizers and other 

the output estimates and the input of symbol decision stage 20 filtering elements operate as "summers" and in accordance 

216. For systems that employ channel coding, an alternative with the well-known Central Limit Theorem cause arbitrary 

way to generate the error signal is to reencode the decoded probability distributions to become Gaussian. The indepen- 

output of channel decoding stage 218 and form an ideal dence assumption also holds for digital cable systems 

modulation signal representing this corrected encoded data. because they generally incorporate large interleavers or error 

Hie error signal is then the difference between this encoded 25 randomizers that break up burst errors (or distortions that 

data and the input of symbol decision stage 216. might cause burst or cause dependence of the errors) and 

FIG. 6A is a block diagram describing various modes of randomize them resulting in independence of error events, 
taking advantage of error signal 217 which represents Bit Error Rate Based on Error Signal Analysis 
ingress interference in communication channel 104. One One embodiment of the present invention makes use of 
possibility is direct display of the error signal in the time 30 the above two assumptions and approximates the signal-to 
domain. An operator able to view this display may be able noise ratio (SNR) to be the cluster variance. Each modula- 
to recognize the source of ingress error as for example, 60 tion scheme has an associated BER vs SNR graph which is 
Hz line noise. Burst errors or other time varying phenomena used to estimate the BER from the measured cluster van- 
may be identified. Alternatively, it may be useful to convert ance. The signal margin may be estimated similarly. For 
the signal into the frequency domain (block 602) using 35 systems that make use of error correction, the post-decoded 
well-known techniques prior to display. This provides the bit error rate (DBER) is estimated by first computing the 
equivalent of a spectrum analyzer on an unoccupied channel. pre-decoded error rate from the cluster variance and apply - 
A frequency domain display may show the presence of ing the input vs output error rate curves for a given decoder 
ingress interference from an AM or FM broadcast signal to compute the output error rate of the decoder. Curves for 
leaking into the cable system, for example. 40 a variety of error correcting codes are efficiently stored in a 

The ingress display can also be normalized to the signal look-up table within block 604. There is no need to simulate 

power to display the power spectral density of the ingress the operation of the decoder. 

signal relative to the signal power. The signal power may be Coarse System Margin Approximation Based on Error Sig- 

obtained at the output of adaptive equalizer 214 or 236. This nal Analysis 

signal power then is used to generate a reference power level 45 Avoiding simulation of decoder operation produces great 

on the frequency display (e.g., OdBc or decibels relative to benefits in determining margin under low error rate condi- 

the signal or carrier power). The ingress power is then tions. This is because measuring low error rates can be very 

plotted (block 604) as a spectral density relative to the signal long and tedious. With powerful error-correction coding, 

power (e.g., -40 dBc/Hz) or as a power in the resolution reasonably short measurements will produce either zero 

bandwidth of the display. For example, if the resolution 50 errors above the SNR threshold or a 50% bit error rate below 

bandwidth is 1 kHz, the -40 dBc/Hz above is multiplied by the SNR threshold, due to the breakdown of the error- 

1000 Hz (30 dB) to obtain the ingress power of -10 dBc. correction coding scheme. This gives very little indication of 

The ingress display can be further corrected for adjust- available margin. However, subtracting a required system 

ments made in the equalization by adjusting in the frequency SNR from the SNR estimate derived from the cluster 

domain for the influence of the equalizer spectral response 55 variance yields a useful estimate of the available margin. A 

on the ingress signals to determine the original ingress technique for obtaining a more accurate estimate of system 

power prior to equalization processing. margin is discussed below. 

Alternatively, statistics of error signal 217 may be com- Precise Evaluation of System Margin by Analysis of Error 

puted (block 604) in relation to the symbol stream output Signal 

from symbol decision stage 216 to develop indicia of line 60 A more refined estimate of system margin may be 

quality. One extremely useful measure of quality is the obtained in accordance with one embodiment of the present 

"cluster variance" (CV). The cluster variance is the ratio of invention. This is done by separating the component of the 

the variance of the power in the error signal to the power in error signal which is caused by broadband noise from the 

the symbol stream. components caused by narrow-band interferers, residual 

Cluster Variance Based on Error Signal Analysis 65 amplitude modulation, and residual phase modulation. 

The cluster variance is measured around the demodulated FIG. 7C depicts a constellation point 756 as affected by 

signal constellation point (e.g., at the output of equalizer 214 broadband noise only. Ideally, constellation point 756 would 
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be a dimensionless point but the effect of broadband noise is 
to diffuse the received constellation point into a disk. FIG. 
7D depicts constellation point 756 as affected by a narrow- 
band interferer. The effect of the narrowband interferer is to 
create a ring 758 having the ideal location of constellation 5 
point 756 as a ceoter. 

FIG. 7E depicts constellation point 756 as affected by 
residual phase modulation. FIG. 7E assumes that constella- 
tion point 758 is in either the upper right or lower left 
quadrant of a constellation such as the constellation in FIG. 10 
2B and located on a line through the origin at roughly a 45 
degree angle to the axes. The effect of residual phase 
modulation is to disperse the received constellation point in 
a direction tangent to a circle including constellation point 
756 and having the origin of FIG. 2B as a center. It should 15 
be noted that the residual phase modulation of FIG. 7E will 
not include errors corrected for by carrier recovery stages 
222 or 238. 

FIG. 7F depicts constellation point 756 as effected by 
residual amplitude modulation. The same constellation 
geometry assumptions made in reference to FIG. 7E also 
apply to FIG. 7F. The effect of residual amplitude modula- 
tion is to disperse the received constellation point in a 
direction perpendicular to a circle including constellation 
point 756 and having the origin of FIG. 2B as a center. It 
should be noted that the residual amplitude modulation of 
FIG. 7F will not include the errors corrected for by AGC 
stage 204. 

For an actual received signal, the errors of FIGS. 7C-7F 
will be combined. The present invention provides a tech- 
nique for isolating these errors. 

FIG. 7G is flowchart describing steps of analyzing an 
error signal to estimate system margin and obtain other 
system information in accordance with one embodiment of 



to FIG. 7J, constellation point 756 lies at a distance r from 
the origin and a vector 778 from origin 776 to constellation 
point 756 is at an angle 6 to a horizontal axis 780. 
Accordingly, the following relationships will hold: 



E{c,€,}*oV(r cos8) 2 o 2 AJV +<r sine) 2 © 2 ^, 
E{iQ*v}=<?Q+{r sine) 2 a 2 AA< +(r cos8) 2 o 2 PW 



(3) 
(2) 
(3) 



20 



where o 2 / and a 2 ^ are the in-phase and quadrature compo- 
nents of the broadband noise respectively, and are each equal 
to cr 2 ;/! The preferred embodiment solves for o 2 ,^ o 2 ^^ 
and cr 2 ^ using the above equations. 

At step 782, the preferred embodiment uses its knowledge 

of 

^ am* pat* an d cr 2 ^ to compute the bit error rate. 
Referring now to FIG. 7K, a decision region 784 surround- 
ing ideal constellation point 756 is shown. Constellation 
point 756 is shown with a combination of mostly residual 
phase modulation with some broadband noise added. As 
long as the constellation point representing the received 
noisy signal falls within the boundaries of decision region 
784, the originally transmitted symbol will be correctly 
identified by the receiver. To evaluate the probability of an 
25 error occurring, the preferred embodiment runs a Monte 
Carlo simulation where the ideal transmitted constellation 
point is summed with pseudo-randomly chosen values cor- 
responding to simulated broadband noise, residual ampli- 
tude modulation, and residual frequency modulation with 
30 distributions defined in part by the o 2 ^, o 2 ^^ and o 2 ^ 
values identified in step 780. The broadband noise and 
residual phase modulations are assumed to have Gaussian 
distributions variances o 2 ^ and cF PM respectively. The 
residual amplitude modulation is modeled as an (e.g., 120 
the present invention. At step 760, the preferred embodiment 35 Hz) periodic random process with an amplitude defined by 
analyzes the spectrum of the received signal as would be ^am- A simulated received point falling outside decision 
present, e.g., on the output of carrier recovery stage 238 or region 784 is considered an error. The bit error rate is 
adaptive equalizer 214 to identify interferers. FIG. 7H determined in this way and displayed. The bit error rate 
depicts such a spectrum 762 showing two narrow band simulation may be refined by taking into account the effects 
interferers 764 and 766. Any narrowband signal having 40 of channel coding by applying simulated channel coding and 
greater than three times the median magnitude of spectrum decoding. 

762 is considered to be an interferer. At step 768, a notch At step 786, the preferred embodiment evaluates the 
filter as shown in FIG. 71 is applied to remove the narrow system margin, determining the amount of extra broadband 
band interference. Referring to FIG. 71, a notch filter noise, residual amplitude modulation, and residual phase 
response 770 is constructed so as to remove the interferers 45 modulation that can be tolerated before the bit error rate rises 



above an unacceptable threshold. TTxis is done by repeating 
the Monte Carlo simulation described in reference to step 
782 for combinations of elevated values of o 2 ^, o 2 ^, and 
cr 2 ^. The combinations for which the bit error rate is 
50 equivalent to the threshold define a three-dimensional sur- 
face which may also be displayed to graphically illustrate to 
the operator how much further system degradation is toler- 
able. 

In summary, the system margin is calculated as the 



identified in step 760. Notch filter response 770 includes 
nulls 772 and 774 corresponding to the frequency and 
bandwidth of narrowband interferers 764 and 766 respec- 
tively. FIG. 71 of course shows ideal filter performance. 

At step 772, the preferred embodiment performs statisti- 
cal analysis on samples from the two-dimensional error 
signal 217 corresponding to a selected constellation point, 
an in-phase error correlation signal E {€,€,} corresponding 
to variance along the in-phase (horizontal) axis, a quadrature 

error correlation signal E {e G €^} corresponding to variance 55 amount of noise power that can be added while maintaining 



along the quadrature (vertical) axis, and an error cross- 
correlation signal E {e^} corresponding to the cross- 
correlation between error samples along the in-phase and 
quadrature axes. 

At step 774, the preferred embodiment determines cr 2 ^, 
corresponding to the amount of residual AM, o 2 ^ corre- 
sponding to the amount of residual PM, and cr 2 ^ correspond- 
ing to the amount of wideband noise based on the correlation 
signals computed in 772. Step 774 takes advantage of a 
relationship among E{e^}, Efe^}, Efe^}, o 2 ^, o 2 ^, 
and o 2 ^ determined by the position of the selected constel- 
lation point in relationship to an origin 776. Referring now 



BER at is threshold value while- considering the cumulative 
effect of other impairments. The present invention, however, 
also provides estimates of the degree to which non-noise 
impairments may be increased before bringing BER above 
60 threshold. 

Demodulation of the error signal (block 606) may also be 
useful. Converting the demodulated signal to audio may 
allow ready identification of broadcast interference. 
Alternatively, the demodulated error signal may be used to 
65 recover one or more deliberately injected in-band ingress 
carriers. If the level of these carriers is sufficiently low (e.g., 
-27 dBc for 64-QAM), the digital carrier is not appreciably 
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disturbed. FIG. 6B depicts communications channel 104 
with an additional ingress signal being injected by a second 
digital transmitter 620. 

Evaluation of System Margin by Applying a Synthetic 
Impairment 

FIGS. 7A-7B are diagrams illustrating bow system mar- 
gin may be tested non-invasively in accordance with on 
embodiment of the present invention. FIG. 7A depicts one 
approach to testing system margin. A summer 702 is inserted 
prior to adaptive equalizer 214 to add in an impairment 
signal. The output of symbol decision stage 216 can then be 
analyzed to measure the deterioration in performance that 
results from the impairment signal in the manner described 
in connection to FIG. 6A. 

An impairment generator 704 generates the impairment 
from error signal 217 generated by symbol decision stage 
216. One approach is simply to amplify error signal 217 
which represents the existing interferes within impairment 
generator 704 until the bit error rate becomes unsatisfactory. 
The degree of amplification necessary to produce this level 
of signal quality degradation then represents how much 
more of the current interference can be tolerated. 
Alternatively, filtering could be applied within impairment 
generator 704 to isolate a particular type of interference to 
be re-applied to the signal and determine how much more of 
this interference could be tolerated. Another alternative is to 
apply additive noise through summer 702 and determine 
how much noise is necessary to disrupt communication. One 
or more reflections may be simulated by using impairment 
generator 704 to delay a reconstruction of the transmitted 
signal based on the symbol estimates. Also, if desired, only 
a particular spectral portion of the reconstruction may be 
delayed in this way. 

FIG. 7B depicts an alternative approach to estimating 
system margin. An impairment processor 750 may apply 
various form of signal impairment to a reconstruction of the 
transmitted signal based on the output of symbol decision 
stage 216. If forward error correction is being used, impair- 
ment processor 750 may instead use the corrected data as 
input and re-apply the relevant code prior to further pro- 
cessing. Impairment processor 750 may apply additive 
impairments as in FIG. 7Aor may apply a simulated channel 
response. The output of impairment processor 750 is applied 
to an equalizer 752 similar to equalizer 214. (The updating 
of the weights of equalizer 752 is not shown.) Equalizer 752 
thus seeks to correct for the impairments introduced by 
impairment processor 750. A symbol decision stage 754 then 
makes symbol decisions based on the output of equalizer 
752. Together equalizer 752 and symbol decision stage 754 
constitute an auxiliary demodulator. The effects of the 
impairments introduced by impairment processor 750 and 
analyzed based on error and symbol estimate outputs of 
symbol decision stage 754 as described in reference to FIG. 
6A. 

Impairment processor 750 may apply additive interfer- 
ence based on error signal 217 in a manner similar to that 
discussed in reference to FIG. 7A. As the interference level 
is varied, the error and symbol estimate outputs of symbol 
decision stage 754 may be monitored for unacceptable 
performance deterioration. 

Impairment processor 750 may also apply a digital filter 
that simulates an intensification of an existing distortion. 
Impairment processor 750 may determine the existing dis- 
tortion by taking the pseudoinverse of the weights 221 
applied by weight updating stage 220. If a particular 
distortion, e.g, a reflection is to be applied, it may be isolated 
by use of simple deconvolution, homomorphic signal 



78,730 Bl 

16 

processing, or parametric signal processing as was discussed 
with reference to FIG. 5A. Varying the intensity level of the 
simulated reflection or reflections while monitoring the 
outputs of symbol decision stage 754 helps determine how 
5 much more impairment can be tolerated without disrupting 
communications. 

Analysis of Receiver Loop Parameters 

FIGS. 8A-8C are diagrams illustrating how parameters 
relating to residual modulation and symbol timing variation 

10 may be measured in accordance with the invention. FIG. 8 A 
depicts how parameters relating to residual amplitude modu- 
lation may be measured and displayed in accordance with 
one embodiment of the present invention. Residual ampli- 
tude modulation may be caused by either transmitter or 

15 receiver components. Referring to receiver 106, residual 
amplitude modulation will make itself apparent in several 
ways. One way is variation in gain control signal 205 which 
will adjust so as to remove residual amplitude modulation 
within its loop bandwidth. The variation in the overall gain 

20 of equalizer 214 as controlled by weight update stage 220 
will operate to remove amplitude variations seen within 
adaptive equalizer 214. Amplitude modulation that is not 
removed by the operation of AGC stage 204 and adaptive 
equalizer 214 will appear as a radial component in error 

25 signal 217 generated by symbol decision stage 216. 

In order to generate information about amplitude 
modulation, an amplitude modulation processing stage 802 
receives gain control signal 205, the weights generated by 
eight updated stage 220 and error signal 217 generated by 

30 symbol decision stage 216. Prior to further processing, gain 
control signal 205 is filtered to remove the known dynamics 
of the gain control loop. The received weights are processed 
to isolate the time varying portion of the overall gain of 
equalizer 214. The cutoff frequency of this isolated gain 

35 response is determined responsive to the loop bandwidth of 
the AGC loop formed by gain stage 202 and automatic gain 
control stage 204. Two-dimensional error signal 217 is 
processed to isolate a radial component. Amplitude modu- 
lation stage 802 then determines, in the time domain, a 

40 weighted sum of gain control signal 205, the time varying 
portion of the gain of equalizer 214, and radial component 
of error signal 217 to obtain a time domain representation of 
the residual amplitude modulation of the received signal 
input to gain stage 202. One of skill in the art will understand 

45 how to adjust the weights appropriate to the parameters of 
adaptive equalizer 214 and symbol decision stage 216. 

Once a time domain representation of the residual ampli- 
tude modulation is generated, it may be displayed directly. 
Alternatively, it may be converted to the frequency domain 

50 using well known techniques prior to display as an ampli- 
tude modulation spectrum. This type of display is highly 
useful in isolating, e.g., hum caused by poor isolation 
between power line noise and transmitter or receiver com- 
ponents. Another useful parameter that may be displayed 

55 numerically in the percentage of amplitude modulation 
which is defined as the amplitude of the time domain 
amplitude modulation signal divided by the overall received 
signal amplitude, which is a fixed parameter of the receiver. 
FIG. 8B depicts how parameters relating to frequency 

60 and/or phase modulation may be determined in accordance 
with one embodiment of the present invention. Referring 
again to receiver 106, residual frequency and/or phase 
modulation will show up in despinning control signal 223 as 
it seeks to correct for the angular component of error signal 

65 217. Equalizer 214 will also seek to correct for residual 
frequency and/or phase modulation and variation through 
variation in the phase component of the weights generated 
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by equalizer weight update stage 220. However, the remain- 
ing uncompensated angular component of error signal 217 
will also represent a component of frequency and/or phase 
modulation. 

In accordance with one embodiment of the present 5 
invention, a frequency modulation processing stage 804 
receives despinning control signal 223, the weights gener- 
ated by weight update stage 220, and error signal 217, Prior 
to further processing, despinning control signal is filtered to 
remove the known characteristics of the carrier recovery 10 
loop. Frequency modulation processing stage 804 isolates 
the time varying phase component of the weights and the 
angular component of error signal 217 and generates a 
weighted sum of these components with despinning control 
signal 223. The weighted sum serves as a time domain 15 
representation of instantaneous phase error. 

This phase error signal can be displayed directly in the 
time domain, enabling the operator to recognize periodicity 
or other patterns. The preferred embodiment also has the 
capability of converting the phase error signal to the fre- 20 
quency domain. The carrier signal, whose amplitude and 
frequency are fixed receiver parameters, can be superim- 
posed over a frequency domain display of the phase error to 
provide convenient read out of phase noise level in dBc/Hz 
(decibels below carrier per Hertz). The phase error signal 25 
can also be differentiated to provide a time domain display 
of frequency error relative to expected carrier frequency. 
Low frequency variations in frequency error represent drift 
of carrier frequency. A stability percentage over a given time 
may also be computed and displayed. An average frequency 30 
error can also be displayed numerically as a fixed offset 
indicating, e.g., transmitter component variation. Note that 
some receiver architectures do not downconvert received 
signals all the way to DC. In these architectures, downcon- 
verted frequency offset should be removed prior to display 35 
of frequency error. A frequency domain representation of the 
differentiated phase error signal will help show sources of 
residual FM. 

FIG. 8C depicts how baud rate and baud rate variation 
may be analyzed and displayed in accordance with one 40 
embodiment of the present invention. Variations in baud 
time will appear within receiver 106 via sampling rate 
control signal 213 which operates to synchronize sampling 
and via the variation in delay through equalizer 24. A baud 
timing processing stage 806 receives sampling rate control 45 
signal 213 and the weights generated by equalizer weight 
update stage 220. Prior to further processing, baud timing 
processing stage 806 filters sampling rate control signal 213 
to remove the known dynamics of the baud timing recovery 
loop. Baud timing processing stage 806 isolates the time 50 
varying delay component of the response of equalizer 214 as 
represented by the weights and obtains a weighted sum of 
sampling rate control signal 213 and the time domain 
representation of the isolated delay component to form a 
timing signal. The so-obtained timing display may be dis- 55 
played numerically as the current baud rate or graphically to 
show variation in baud rate over time. Baud rate variations 
would be caused by errors or problems in the transmitter. A 
baud rate stability percentage over a given time may also be 
displayed. 60 

Of course, in receiver architectures that do not include an 
adaptive equalizer, estimates of the above-described modu- 
lation and timing parameters could be obtained from the 
various control loop parameters alone. Equalizer weights 
and error signals would be considered. 65 

In a preferred embodiment, the invention is sufficiently 
modular and programmable to allow customization of the 
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demodulator, system parameters, modulation types, displays 
and diagnostics to meet specific system requirements. 

In another preferred embodiment, the invention is pow- 
ered by a battery, increasing the portability of the system. 

Variations on the methods and apparatus disclosed above 
will be readily apparent to those of ordinary skill in the art. 
Accordingly, particulars disclosed above are intended not as 
limitations of the invention. 
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We claim: 

1. In a digital communication system wherein digital data 
to be communicated over a channel is represented by 
sequential selection of modulation components represent- 
able as points of a two-dimensional constellation wherein a 
horizontal position of a given point indicates a magnitude of 
an in-phase component of a modulated signal transmitted via 
said channel and a vertical position of said given point 
indicates a magnitude of a quadrature component of said 
modulated signal, a method of diagnosing said channel 
comprising the steps of: 

a) receiving said modulated signal as corrupted by said 
channel; 

b) estimating an error signal between said modulated 
signal as transmitted and said received modulated sig- 
nal; 

c) estimating, using said error signal, an in-phase error 
correlation signal, a quadrature error correlation signal, 
and a cross-correlation error signal; and 

d) estimating a component of said error signal responsive 
to said in-phase error correlation signal, said quadrature 
error correlation signal, and said cross-correlation error 
signal, and a position of said given constellation point 
within said constellation. 

2. The method of claim 1 wherein step d) further com- . 
prises estimating a residual phase modulation component of 
said error signal responsive to said in-phase error correlation 
signal, said quadrature error correlation signal, and said 
cross-correlation error signal, and the position of said given 
constellation point within said constellation. 

3. The method of claim 1 wherein step d) further com- 
prises estimating a residual amplitude modulation compo- 
nent of said error signal responsive to said in-phase error 
correlation signal, said quadrature error correlation signal, 
and said cross-correlation error signal, and the position of 
said given constellation point within said constellation. 

4. The method of claim 1 wherein step d) further com- 
prises estimating an additive noise component of said error 
signal responsive to said in-phase error correlation signal, 



07/13/2004, EAST Version: 1.4.1 



US 6,278, 

19 

said quadrature error correlation signal, and said cross- 
correlation error signal, and the position of said given 
constellation point within said constellation. 

5. The method of claim 1 further comprising the steps of: 
registering a spectrum of said error signal; 5 
identifying a narrowband interfering signal from said 

spectrum; and 

prior to said c) step, applying a notch filter to said error 
signal to remove the narrowband interfering signal 
from said error signal. 

6. The method of claim 1 wherein step d) further com- 
prises estimating a plurality of components of said error 
signal responsive to said in-phase error correlation signal, 
said quadrature error correlation signal, and said cross- 
correlation error signal, and the position of said given 
constellation point within said constellation, said plurality of 
components including a residual phase modulation 
component, a residual amplitude modulation component and 

an additive noise component. 2Q 

7. The method of claim 6 further comprising the step of 
estimating a bit error rate response to said additive noise 
component, said amplitude modulation component, and said 
residual phase modulation component. 

8. The method of claim 6 further comprising the step of 25 
estimating tolerable additional levels of additive noise, 
residual amplitude modulation, and residual phase modula- 
tion. 

9. The method of claim 7 wherein said bit error rate 
estimating step further comprises performing a Monte Carlo 3Q 
simulation of said additive noise component, said amplitude 
modulation component, and said residual phase modulation 
component. 

10. The method of claim 8 wherein said tolerable addition 
levels estimating step further comprises performing a Monte 3S 
Carlo simulation of said additive noise component, said 
amplitude modulation component, and said residual phase 
modulation component. 

11. A computer program product that controls a processor 

to diagnose a channel in a digital communication system in ^ 
which digital data is communicated over said channel as 
sequential selection of modulation components represent- 
able as points of a two-dimensional constellation wherein a 
horizontal position of a given point indicates a magnitude of 
an in-phase component of a modulated signal transmitted via 45 
said channel and a vertical position of said given point 
indicates a magnitude of a quadrature component of said 
modulated signal, the computer program product compris- 
ing: 

code for estimating an error signal between a received 50 
modulated signal, said received error signal corrupted 
by said channel, and said modulated signal as trans- 
mitted; 

code for estimating, using said error signal, an in-phase 
error correlation signal, a quadrature error correlation 55 
signal, and a cross-correlation error signal; and 

code for estimating a component of said error signal 
responsive to said in-phase error correlation signal, said 
quadrature error correlation signal, and said cross- 
correlation error signal, and a position of said given 60 
constellation point within said constellation. 

12. The computer program product of claim 11 further 
comprising code for estimating an additive noise component 
of said error signal responsive to said in-phase error corre- 
lation signal, said quadrature error correlation signal, and 65 
said cross-correlation error signal, and the position of said 
given constellation point within said constellation. 
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13. The computer program product of claim U further 
comprising code for: 

registering a spectrum of said error signal; 
identifying a narrowband interfering signal from said 
spectrum; and 

prior to said c) step, applying a notch filter to said error 
signal to remove the narrowband interfering signal 
from said error signal. 

14. The computer program product of claim 11 further 
comprising code for estimating a plurality of components of 
said error signal responsive to said in-phase error correlation 
signal, said quadrature error correlation signal, and said 
cross-correlation error signal, and the position of said given 
constellation point within said constellation, said plurality of 
components including a residual phase modulation 
component, a residual amplitude modulation component and 
an additive noise component 

15. The computer program product of claim 14 further 
comprising code for estimating a bit error rate response to 
said additive noise component, said amplitude modulation 
component, and said residual modulation component 

16. The computer program product of claim 14 further 
comprising code for estimating tolerable additional levels of 
additive noise, residual amplitude modulation, and residual 
phase modulation. 

17. The computer program product of claim 15 further 
comprising code for estimating the bit error rate by per- 
forming a Monte Carlo simulation of said noise component, 
said amplitude modulation component, and said residual 
phase modulation component. 

18. The computer program product of claim 15 further 
comprising code for estimating tolerable addition levels of 
additive noise, residual amplitude modulation, and residual 
phase modulation by performing a Monte Carlo simulation 
of said noise component, said amplitude modulation 
component, and said residual phase modulation component. 

19. An arrangement for diagnosing a channel for use in 
connection with a digital communication system for com- 
municating digital data to be communicated over said chan- 
nel as sequential selection of modulation components rep- 
resentable as points of a two-dimensional constellation 
wherein a horizontal position of a given point indicates a 
magnitude of an in-phase component of a modulated signal 
transmitted via said channel and a vertical position of said 
point indicates a magnitude of a quadrature component of 
said modulated signal, said arrangement comprising: 

a sampling stage configured to receive the modulated 
signal and generate a digital signal synchronized to a 
plurality of symbol boundaries in said received modu- 
lated signal; 

an equalizer configured to received the digital signal and 
generate an equalized modulated signal based on the 
digital signal; 

a symbol decision stage configured to received the equal- 
ized modulated signal and generate a plurality of con- 
stellation points corresponding thereto, said symbol 
decision stage further configured to estimate an error 
signal between the received modulated signal and the 
modulated signal as transmitted; and 
a processor configured to 

estimate, using said error signal, an in-phase error 
correlation signal, a quadrature error correlation 
signal, and a cross-correlation error signal; and 
estimate a component of said error signal responsive to 
said in-phase error correlation signal, said quadrature 
error correlation signal, and said cross-correlation 
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error signal, and a position of said given constella- 
tion point within said constellation. 

20. The arrangement of claim 19 wherein the sampling 
stage is operably connected to receive the modulated signal 
from a tuning/downconversion circuit in analog format. 

21. The arrangement of claim 19 wherein the sampling 
stage is a digital resampler operably connected to receive the 
modulated signal from a quadrature direct conversion stage 
in digitally sampled format. 

22. The arrangement of claim 19 wherein the equalizer 
comprises an adaptive equalizer configured to correct for 
distortion in the channel based on the error signal. 

23. In a digital communication system wherein digital 
data to be communicated over a channel is represented by 
sequential selection of modulation components represent- 
able as points of a two-dimensional constellation wherein a 
horizontal position of a given point indicates a magnitude of 
an in-phase component of a modulated signal transmitted via 
said channel and a vertical position of said given point 
indicates a magnitude of a quadrature component of said 
modulated signal, a method of diagnosing said channel 
comprising the steps of: 

a) receiving said modulated signal as corrupted by said 
channel; 
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b) estimating an error signal between said modulated 
signal as transmitted and said received modulated sig- 
nal; 

c) estimating a cluster variance responsive to said error 
signal; and 

d) estimating a bit error rate from said cluster variance. 

24. The method of claim 23 wherein step c) further 
comprises: 

estimating a power of said received modulated signal; 
from said error signal, estimating a sum of a power of an 

in-phase component of said error signal and a power of 

a quadrature component of said error signal; 
estimating a cluster variance as a ratio of said power of 

said received modulated signal and said power of said 

sum. 

25. The method of claim 23 wherein said bit error rate is 
a pre-decoded error rate. 

26. The method of claim 25 further comprising the step of 
converting said pre-decoded bit error rate to a post-decoded 
bit error rate according to an error correction coding scheme 
employed by said channel. 
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